
Mesh.MshFileVersion = 1;
Mesh.MshFileVersion = 1;
a=-1;
b=1;
c=-1;
d=1;
e=-1;
f=1;
h=0.1;
Point(1) = {a,c,e,h};
Point(2) = {b,c,e,h};
Point(3) = {b,d,e,h};
Point(4) = {a,d,e,h};
Point(5) = {a,c,f,h};
Point(6) = {b,c,f,h};
Point(7) = {b,d,f,h};
Point(8) = {a,d,f,h};

Line(1) = {1,4};
Line(2) = {4,3};
Line(3) = {3,2};
Line(4) = {2,1};
Line(5) = {1,5};
Line(6) = {5,6};
Line(7) = {6,7};
Line(8) = {7,8};
Line(9) = {8,5};
Line(10) = {2,6};
Line(11) = {3,7};
Line(12) = {4,8};
Line Loop(5) = {3,4,1,2};
Line Loop(6) = {-8,-11,-2,12};
Line Loop(7) = {12,9,-5,1};
//Line Loop(8) = {8,12,2,11};
//Line Loop(9) = {7,10,3,11};
//Line Loop(10) = {6,10,4,5};
Plane Surface(6) = {5};
Plane Surface(7) = {6};
Plane Surface(8) = {7};
Physical Line(1) = {1};
Physical Line(2) = {2};
Physical Line(3) = {3};
Physical Line(4) = {4};
Physical Line(5) = {5};
Physical Line(6) = {6};
Physical Line(7) = {7};
Physical Line(8) = {8};
Physical Line(9) = {9};
Physical Line(10) = {10};
Physical Line(11) = {11};
Physical Line(12) = {12};
Physical Surface(6) = {6};
Physical Surface(7) = {7};
Physical Surface(8) = {8};
nx = 1/h;
ny = 1/h;
nz = 1/h;
Transfinite Line {1,3} = ny + 1 Using Progression 1.0;
Transfinite Line {2,4} = nx + 1 Using Progression 1.0;
Transfinite Line {4,6} = nz + 1 Using Progression 1.0;

Transfinite Surface {6} = {1,2,3,4};
Transfinite Surface {7} = {5,6,7,8};
Transfinite Surface {8} = {1,5,8,4};
Recombine Surface {6};
Recombine Surface {7};
Recombine Surface {8};
Line Loop(13) = {11,-7,-10,-3};
Ruled Surface(14) = {13};
Line Loop(15) = {6,-10,4,5};
Ruled Surface(16) = {15};
Line Loop(17) = {9,6,7,8};
Ruled Surface(18) = {17};
